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Abstract 

We study the dissipative quantum Duffing oscillator in the deep quantum 
regime with two different approaches: The ffist is based on the exact Floquet 
states of the linear oscillator and the nonlinearity is treated perturbatively. 
It well describes the nonlinear oscillator dynamics away from resonance. The 
second, in contrast, is applicable at and in the vicinity of a A^-photon reso- 
nance and it exploits quasi-degenerate perturbation theory for the nonlinear 
oscillator in Floquet space. It is perturbative both in driving and nonlin- 
earity. A combination of both approaches yields the possibility to cover the 
whole range of driving frequencies. As an example we discuss the dissipative 
dynamics of the Duffing oscillator near and at the one-photon resonance. 

Keywords: Nonlinear quantum oscillator, Floquet theory 
1. Introduction 

Classical nonlinear systems show interesting phenomena like bistability, 
frequency doubling and nonlinear response and cover a wide range of appli- 
cability and various physical realizations H]- In the last years it has been 
possible to build nonlinear devices which can potentially reach the quan- 
tum regime. These are for example cavities incorporating a Josephson junc- 
tion I3L4L SQUIDs used as bifurcation amplifiers to improve qubit read-out 
[sl, 0,^7, ij or nanomechanical resonators jo], |lo[. Recently, a novel class 
of devices combining SQUIDS and resonators has been demonstrated. For 
example, sensitive detection of the position of a micromechanical resonator 
embedded in a nonlinear, strongly damped DC-SQUID has been achieved 



In the deep quantum regime, where bistability is no longer observed. 
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there has been to date no experimental operation of nonhnear driven oscil- 
lators to our knowledge. 

From the theoretical side, semiclassical approaches have been used to describe 
the situation where the underlying classical bistability still plays a dominant 
role. A DC-SQUID embedded in a cavity allowing displacement detection 
and cooling was analyzed in lif. Composed qubit- Josephson bifurcation 



amplifier systems have been considered in Il3l . |14| . Dynamical tunneling in 
a Duffing oscillator was accounted for in [l5| within a semiclassical WKB 
scheme, while in 16|, ll7|, |l8| a Wigner function analysis near the bifurcation 
point is put forward. 

The behaviour of the Duffing oscillator (DO) in the deep quantum regime 



has attracted recently lot of interest. In particular Rigo et al. [19| demon- 
strated, based on a quantum diffusion model, that in the steady state the 
quantum DO does not exhibit any bistability or hysteresis. It was also shown 
that the response of the Duffing oscillator displays antiresonant dips and res- 
onant peaks [20|, l2l|, |22|, |23| depending on the frequency of the driving field, 
originating from special degeneracies of the eigenenergy spectrum of the non- 
linear oscillator [20']. While the antiresonances persist in the presence of a 
weak Ohmic bath, for strong damping the nonlinear response turns to a res- 
onant behaviour, namely the one of a linear oscillator at a shifted frequency 



22(1 . Finally, recently Nakano et al. [2J] looked at the composed qubit-DO 



dynamics during read-out process. 

In this work we investigate the deep quantum limit of the quantum Duffing 
oscillator and present two different approaches covering different parameter 
regimes. The first approach is based on the exact Floquet energies and states 
of the driven linear oscillator with the nonlinearity treated perturbatively. 
As there is no restriction on the driving amplitude, this scheme can also be 
applied to the regime where the driving amplitude is larger than the non- 
linearity. The second approach treats both the driving and the nonlinearity 
perturbatively. It is applicable for driving frequencies which can resonantly 
excite two states of the nonlinear oscillator, requiring that the driving cannot 
overcome the nonlinearity. In general a combination of both approaches al- 
lows to cover the whole range of driving frequencies. Exemplarily we consider 
the dynamics of the Duffing oscillator near the one-photon resonance, where 
the oscillator dynamics is described analytically. As in 20|, |2l|, |22|, |23 



we 



obtain that for weak dissipation the amplitude of the oscillations displays an 
antiresonance rather than a resonance. We find a characteristic asymmetry 
of the antiresonance lineshape. In contrast to 20|, |2l|, |22|, |23|, our analytic 
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results are obtained without applying a rotating wave approximation (RWA) 
on the Duffing oscillator. 

The paper is organized as follows: In section |2] we introduce the Hamil- 
tonian of the non-dissipative Duffing oscillator and the two Floquet based 
approximation schemes to treat it. The energy spectrum and eigenstates of 
the non-dissipative system are calculated with the two different schemes in 
section [3] and section [Hand both approaches are compared in section O After- 
wards dissipative effects are included within a Born-Markov-Floquet master 
equation in section |6l Section [7] addresses the special case of the one-photon 
resonance including dissipative effects. In section [8] conclusions are drawn. 



2. Quantum DufRng oscillator 

A quantum Duffing oscillator is described by the Hamiltonian: 



p2 

y 



(1) 



2M ■ 2 " ■ 4 

where M and Q are the mass and frequency of the Duffing oscillator which 
is driven by a monochromatic field of amplitude F and frequency Uex- For 

later convenience we introduce the oscillator length yo : 



h 



In the fol- 



lowing we will consider the case of hard nonlinearities, a > 0, such that the 
undriven potential is monostable. 

To treat the quantum Duffing oscillator problem we observe that the Hamil- 
tonian can be rewritten as: 



a 



HLoit) + jy' 
Htiho + yFcos{uext), 



(2a) 
(2b) 



where Hi,o{t) describes a driven linear oscillator, while f^NLO is the Hamilto- 
nian of an undriven nonlinear oscillator. Due to the periodic driving Floquet 
theory [25|, [iol, 27, 28|, reviewed in Appendix A can be applied. In particu- 
lar Floquet theorem states that solutions of the time-dependent Schrodinger 
equation for Hj)o(t) are of the form 

lV'iW)=exp(-^e,t/n)|0,(t)), (3) 

where \(pj(t)) = \(pj(t + T^^^)). The quasienergies ej and Floquet states 
|<^j(t)), respectively, are eigenvalues and eigenf unctions of the Floquet Hamil- 
tonian 'HBoit) = HBoit) — ihdt- As discussed in Appendix A |0j,n(i)) = 



3 



\(f)j{t)) exp{—inujext) yields a physically equivalent solution but with shifted 
quasienergy ej^n = — nhuex- Eqs. (12a|l and pbj) suggest two different ap- 
proaches, shown in Figure [H to solve the eigenvalue problem described by 
Eqs. (IA.13|) and flA.14|) : In the first one, called App I, starting point are the 



Driving 

App I I 

Nonlinearity i Nonlinearity 

App II \ 

-f^NLO --f^DO {t) 

Driving 

Figure 1: Different procedures to incorporate driving and nonlinearity. In App I starting 
point are the exact Floquet states and eigenenergies of the driven linear oscillator HLo{t)- 
The nonlinearity is the perturbation. In App II the driving is a perturbation expressed on 
the basis of the Floquet states of the undriven nonlinear oscillator Hnlo- 

exact Floquet states and eigenenergies of the driven linear oscillator Hi^oit), 
see Eq. ( |2al) . The nonlinearity is treated as a perturbation. A similar prob- 



lem was considered by Tittonen et al. 29|. This approach is convenient if 
the Floquet states of the time-dependent Hamiltonian are known. 
For the driven harmonic oscillator they have been derived by Husimi and 
Perelomov 



30, |3l| and are given in [Appendix B 



In the second approach, which we call App II, one considers as unperturbed 
system the undriven nonlinear oscillator (NLO) and the driving is the per- 
turbation, see Eq. ( l2b|) . 

As we shall see, the different ways of treating the infinite dimensional Floquet 
Hamiltonian result in crucial differences when evaluating observables of the 
Duffing oscillator. 



3. Perturbation theory for a time-periodic Hamiltonian with time- 
independent perturbation 

The starting point of the perturbative treatment App I is the Floquet 
equation for the full Floquet Hamiltonian l-L in the extended Hilbert space 
7^® r, see Eq. f[03|) . 

= ej^rn\(t)j,m)), (4) 
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where l-i = i-Lo+Va- Moreover, the Floquet states of the Floquet Hamiltonian 
iio satisfying the eigenvalue equation (1A.13P are known, see e.g. Eq. (IB. lap 
and Eq. flRi]) : 



^O|0i,m))o = (^fj<Pj,m))o- (5) 

We look for an expression of ej^m and \(pj,m)) in first order in Va- Hence we 
introduce the first order corrections e^-^^ and \(j)j,m))i as: 



|0j,m)) = |0j,m))o + |0j,m))l- 

Because the perturbation is time- independent, it is diagonal in the Hilbert 
space T. AdditionalUy, we introduce the Fourier coefficients: 

As in the case of conventional stationary perturbation theory, the perturbed 
states are written as a linear combination of the unperturbed states: 

\(l>,,m)) = |0,,J)O+ E d'^."))0' (8) 

where (i, n) denotes the couple of quantum numbers i and n. Inserting ansatz 
Eq. ^ in the Floquet equation (jlj) we obtain: 

(?io + V;-65.°i-e5.^i)(|0,,^))o + |0,,J)i) = 0. (9) 

Because e^^^ and |0j,m))o solve the Floquet equation for T-Lq the last equation 
reduces to: 

t^a-e5.^i)|0,,^))o+(Ho-65°i)|0,,„))i = 0. (10) 



(1) 



This equation allows to determine the modification to the quasienergy 
and the actual form of the coefficients c^™. To calculate ef^^ we multiply it 
from the left with o((0i,m|- It follows: 



o{{4>j,m\Va\(pj,m))o- (H) 



To determine the coefficients for the states we multiply Eq. fllUp from the left 
with o((0fc,n|, where the couple {k, n) ^ (j, m). Moreover we exclude the case 
of degenerate quasienergies and impose: e"^\ ^ ef^. In case of e[,°^ = e^°^ 
degenerate perturbation theory should be applied. We obtain: 



= o((0MiK.i0.,jo+(4°i-4!l)cr' 



(12) 



yielding: 



(n—rri) 



(13) 



If we set the driving to zero the quasienergy and the states reduce to the 
ones obtained by applying conventional stationary perturbation theory on 
the unforced system. 



3.1. Application to the quantum Duffing oscillator 

We can now determine the actual form of the quasienergy spectrum and 
the corresponding expansion coefficients for the case of the quantum Duffing 
oscillator using as perturbation Va — """'"^ 
defined as: 



y^. The matrix elements Eq. ([7]) 



V 



(n— m) 
kj 



x. 



rftexp(i(n - m)ujext)o{(t)k{t)\Va\(t)j{t))o, (14) 



are given in [Appendix C[ The quasienergies of the quantum Duffing oscil- 



lator, exact in all orders of the driving strength and up to first order in the 
nonlinearity, read: 



a 

AMiul - + 4 



(15) 
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In the limit of no driving Eqs. (jS]) and f llSp yield: 

1\ 3 



lim e,- = hQ j 



F-s>0 



+ o"?/0 U'(J + 1) + 



1 



lim|0,(t)) = lim(t|0,)) 



|j)o + 



(16a) 
(16b) 



4 



2hQ 



\J - 2)o 



, V(j + l)(j + 2)(j + i) 



^iv/.(.-l)C^2)(,-3)^^_^^^^ 



|v/(j + l)(j + 2)(j + 3)(j+4) 



\J + 4)c 



such that the modifications due to the nonlinearity are exactly those obtained 
by conventional stationary perturbation theory 32|, where {| j)o} are the 
eigenstates of the undriven harmonic oscillator. Expanding up to second 
order in the driving amplitude we obtain from Eq. for \4>j(t)) = (t\(j)j)) 
the result: 



a 



-2hn 



l0.+2(t))o(17) 



[y'oU - 1) + lyl^l]^ 



2hn 



v/(j + l)(j + 2)(j + 3)(j+4) 



fVj(j-l)(j-2)(j-3) 



I0i+4(t))o 



l0,-4(t))o 



+ -(2j + l)i/oM| 



+^(j + i)VjTT^^2/o' 



exp ( —i2uJext) exp {i2oJext) 



4 



l0.W)o 
exp(-ia;ex^) _ exp(ia;ex^) 



l0i+iW)o 
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l0.--iW)c 



+ V(j + 3)(j + 2)(j 

+ Vj(J-1)(J-2)^1/o' 
+ ^2/o^cV(j + l)(j + 2) 



exp(-za;e^t) _ exp(+za;ext) 

l_ + 3/lfi —fVjJex + 3/lfi 

exp(-22a;ex'^) _^ exp(z2a;ex^) "' 

h2Uex — —h2Uex — 2MI 



2^/2 

4 -^0^? 



3 / ■/ ■ _ IN rexp(--z2we^ exp(z2a;ext) 

4^° ? V JU J ^ ^2a;,^. + 2h^l -h2cUex + 2hn 



I0i+3(t))o 
l0.-3(t))o 
I0j+2W)O 



where we used the abbreviation = 
states of the hnear oscillator. 



l0.--2(t))o ] , 

M(uj£-n^) \^ji't))o are the Floquet 



4. Perturbative approach for the one-photon resonance 

When the nonlinearity becomes a relevant perturbation to the equidistant 
spectrum of the linear oscillator, it becomes preferable to use the second 
approximation scheme, App II, based on the decomposition Eq. fl2bp . 
In this case it is convenient to express the Floquet Hamiltonian Hbo in the 
composite Hilbert space 71®T spanned by the vectors \i,n)) = \j) ® |n), 
where \j) is an eigenstate of the nonlinear oscillator i^NLO given in Eq. fflGbl) . 
Hence, in this basis the Floquet Hamiltonian of the nonlinear oscillator Hnlo, 
see Eq. (IT^ below at vanishing driving amphtude, is diagonal. In contrast, 
the perturbation Vp = yF cos{uext) is time-dependent and thus non-diagonal 
also in the Hilbert space T. From the relation: 

(O^I'^Dol^,'^')) = {HjiofjT'^^ - ^exndjk^nn', (18) 

where {H-Do)^jl " ^ are the Fourier coefficients of the matrix {j\HDo{t)\k) , it 
follows 

F 

), (19) 

with Ej^n = Ej — hnujex and Ej the energies of the nonlinear oscillator Eq. 
fll6al) . From Eq. (fT9|) is is thus apparent that two eigenstates |j, n)), \k,m)) 
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of Hnlo become degenerate when Ej^n = E^y, i.e. for a driving frequency 
Uex satisfying 

huexin'-n) = Ek-Ey (20) 

Setting N = n' — n one speaks of a A^-photon resonance. From Eq. fll6al) for 
the energies Ej it follows, with k = j + N, 



hUexN = Ej+N -Ej=N 



hn + lay',{N + l + 2j) 



(21) 



In the following we restrict to the one-photon resonance = 1, i.e., the 

quasienergies Ej^n and Ej^i^n+i are degenerate if Uex = + ^"^"4^^"^^ = ^j- 
Moreover, due to the arbitrariness in the choice of the Brillouin zone index 
n, we fix it in the following to the zeroth Brillouin zone, i.e., n = 0. For 
our perturbative treatment we further require that the nonlinearity is large 
enough that if Ej^ is resonant with Ej+i^+i the remaining quasi-energy lev- 
els are off resonance and not involved in the doublet spanned by the two 
degenerate levels. Having this in mind, we have to restrict ourselves to a 
certain range of possible values of Uex, namely to the resonance region such 
that the chosen doublet remains degenerate or almost degenerate, i.e. for the 
one-photon resonance: ~ ^j\ < ^"^^/o- This results from the fact that if 
Ej^Q = Ej^i^^i, the closest lying levels Ej+2,2 and E'j-i^-i are by ^a?/Q away. 
Because of the manifold (doublet) structure of the quasi-energy spectrum. 



we apply in the following Van Vleck perturbation theory [33|, |3J] and treat 
the driving as a small perturbation, i.e. -C J^a?/o ^ Consequently, 
a consistent treatment in App II requires that either contributions are 
neglected if we consider the nonlinearity only up to first order, or that both 
driving and nonlinearity are treated up to second order. As the second order 
in both parameters is very involved, we restrict to the first order in the non- 
linearity and neglect quadratic contributions in the driving strength, as long 
as their reliability cannot be verified within a different approach, i.e. App I. 
Within Van Vleck perturbation theory we construct an effective Floquet 
Hamiltonian = exp(i5')HDO exp{—iS) having the same eigenvalues as 
the original Hamiltonian Hbo and not containing matrix elements connect- 
ing states belonging to different manifolds. Therefore it is block-diagonal 
with all quasi-degenerate energy states in one common block. To determine 
the transformation S and the effective Hamiltonian H^s we write both as a 



9 



power series in the driving: 

S = 5(0) +5^ + 5(2) + ... 
rtcS - ties + ^eff + ^eff + 



(22) 
(23) 



In Appendix D the general formulas for the energies and the states up to 



second order are provided [33|, l3J, l35|, [36 



The zeroth order energies are Ej^ and E'j+i^+i and the corresponding (quasi)- 
degenerate Floquet states are: \j, 0)) and \j + 1, +1)). 

The quasi-degenerate block of the effective Hamiltonian in this basis up to 
second order in the driving strength acquires the form: 



K 



eff 



(1) 



(2) — 



E. 



(1) 



E) 



E. 



+ E. 



(2)++ 



where 



E 



(1) 



((j, 0\Vf\j + 1, +1)) = {{j + 1, +1\Vf\j, 0)) = 



ypE 

2^2' 



and 




Ef>- 
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^(2)++ 
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with 





+ 



njij-l) 
Ej,o ~ -Sj-i.-i Ej^Q — Ej_i^^i 



+ 



Ej,o 



Ej+i^+l - Ej+2,+2 



+ 



Ej+i~i 



(24) 



(25) 



(26) 



E, 



-1 — Ej 



i+2,0 



^i(j) 



8hn 



a(j + 1) 



Ej+i,+i - -Ej,2. 

(27) 



Notice that the unperturbed quasienergies Ej^ and are correct up to 

first order in the nonhnearity. For consistency also n'^{j) has to be treated 
up to first order in a only. 

As shown by Eq. fj28|) below it is essential to determine the eigenenergies of 
T^eff up to second order in E. They are also the eigenenergies of ?^do and 
read: 



^7 



^ [e,,o + + Ef- + ++) ± 1 [(E,,o - i?,+i,+i)' (28) 



+2{Ejfi - Ej+i^+i){E- 



(2)- 



E. 



(2)H 



4E. 



(If 



nl/2 
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The convention ej is chosen such that e„- < e^^ for Uex < , whereas it 
jumps at resonance such that ej > for Uex > ^j- Because the first order 
correction in the driving Ej^'^ enters Eq. fl251) quadratically, a calculation of 
the quasienergies up to first order in F merely yields (when EjQ ^ Ej+i^i) 
the zeroth order results. Consequently to be consistent one has to take into 

^2) (2)++ 

account also the second order corrections Ej and Ej to the energies. 
The eigenstates of the block Eq. fl2^ are determined by: 



|--,o))eff := -sin||j + l,+l)) + cos||j,0)), (29) 
|+,-i))eff := sin^|j,0)) + cos||j + l,+l)), 

where 

tanr^j = rrr . . . (30) 

-(E,,o - + -- - ++) 

In conventional Van Vleck perturbation theory the eigenstates of Hdo ^^i"© 
obtained by applying a back transformation: 

|=Fj,n)) = exp{-iS)\Tj,n))eS- (31) 
Expanding the exponential up to first order we obtain for the eigenstates: 

l"Fi,n)) = |=Fi,n))cfT — '^'S'*-^-' |=Fj,n))cfr, 

= \Tj,n))cS + RyF\Tj,n))cS- (32) 

For reasons given in section [S] and the chosen parameter regime of App II we 
do not determine the second order correction for the states coming from the 
second order contribution to S. 

In the second line of Eq. (!32l) we used the fact that we can express the trans- 
formation S by introducing the reduced resolvent R, allowing a nice connec- 



tion to conventional degenerate perturbation theory as shown in Appendix D 
From Eq. it follows: 

|--o)) = -sin|(l + i?VV)|j + l,+l)) + cos|(l + i?f^)|j,0)), (33) 
= sin|(l + W^)|j,0)) + cos|(l + /2f^)|j + l,+l)), 
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where 

RVf\j,0)) = E -E ^^^^ 

(fc,n)^{(j,0),(i+l,+l)} '"'^ 

|,-l.+l)>+ /"''^;,°> |, + 3.+l)) 

|j-3,+l)) + -^^^%^|j + 3,-l)) 







— 


1 


-3, 


a) 


Ejfi - 




-3,+l 


1 ^3(j 


-3 


a) 



|j-3,-l)) , 



\k,n)){{k,n\VF\j + 



(fc,n)^{(j,0),(j+l,+l)} ^ + ^'+^ 

721(7 + 1) , „ n^f? — 2, a) , 



— Ej+2,0 ~ -^i-2,0 

+^r^^^^ - 2, +2)) + 1^ ' |j + 4, 0)) 

7^3(7 + 1, o) , . , , 
+ l^|j+4,+2)) 



and 



o / 

^3(J,«) = Y^V(j + 3)(j + 2)(j + l). 

The effect of the transformation is to yield a contribution from states outside 
the manifold. Notice that in order to obtain the states to first order in E the 
trigonometric functions sin y and cos y should be expanded in powers of F. 
We conclude this section by mentioning that eigenenergies and eigenstates 
of the Duffing oscillator have been calculated near and at resonance also 



by Peano et al. [22|. However in |22| the nonlinear undriven Hamiltonian 



^NLO = Tm + + fy^ is approximated by H^lo - ^^J + iaVojO + 1), 
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where j is the occupation number operator of the undriven hnear oscillator. 
This approximated Hamiltonian is diagonal in the linear oscillator basis and 
yields the result Eq. fll6ap for the energies of -f^NLO- However, further correc- 
tions of order a contained in the eigenstates fll6bl) are neglected. The results 
of 22| at finite driving can be retained from Eqs. f l28|) and f l33|) by treating 
the driving up to first order, by replacing ni{j) by a/j~+T and by setting 
risU) = 0. 



5. Comparison of the outcomes of the two approaches 

The approximation scheme in Sec. [HI App I, is valid when the quasienergy 
spectrum of the linear oscillator is non-degenerate, i.e., away of a A^-photon 
resonance. In contrast, the perturbative approach of Sec. HI denoted as App 
II, works at best near a A^-photon resonance in the quasienergy spectrum of 
the undriven nonlinear oscillator. Thus a comparison of the outcomes of the 
two approaches is possible in the frequency regime near resonance, i.e. within 



< \ujex ~ ^j\ < ^c^Z/o- AdditionalUy, as the Van Vleck-based approach is 
perturbative in the driving, remember <^ JE'^yo^ ^ comparison requires 
an expansion in F of the results from App I. 

This section is organized as follows: First the energies and then the matrix 
elements of the position operator are compared. 



5.1. Comparison of the quasienergies 
We start with the off resonant case \E 



E. 



h\uJe 



n,\>E^ 



(1) 



and expand the result in Eq. fl28|) up to second order in the driving amplitude 
F: 



E,,o + Ef- +E]'>'/{E 



(1)2 



E. 



(35) 



E. 



(2)++ 



Expanding further for consistency the eigenvalues up to first order in the 
nonlinearity we obtain: 



2 172 r 



+0{F\a' 



2Q 



+ 



(2j + 1) 



fi2 



(36) 



E. 



+ 



2n 



+ 



?,ayt (2j + 3)^]^ 



+ 0{F\a' 
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Inserting yo = ^Jfij (Mil), these are exactly the results obtained from App 
I for e^ o and Cj+i i upon expanding Eq. (fT5i) up to second order in the 
driving amplitude. Consequently, as the quasienergies f|T5|) coincide with the 
quasienergies from App II away from the resonance, we conclude that Eq. 
fl28|) describes the frequency dependence of the quasienergy up to 
over the whole parameter regime, i.e., near and far from resonance. Moreover, 
because the contribution of order to the quasienergies obtained in App 

I is state- independent, see Eq. (1T5|) . it drops when differences of quasienergies 
are considered. In other words the difference of quasienergies coincides in 
both approaches. 

5.2. Matrix element yik{t) 

Due to the agreement (in second order in F) for the quasienergies and the 



disagreement for the Floquet states shown in Appendix E , the question arises 



whether expectation values of observables also differ in the two approaches. 
We shall answer this question in the following at the level of the expectation 
value of the position operator y. 



5.2.1. yikit) in App I 

For the linear oscillator the exact result holds: 



o{Ut)\y\Mt))o 

dy'dy" o {(pi (t) \y'){y'\y\ y") {y" \Mt))o 

dyUy - m)y'4>k{y - m) 



(37) 



1/0 



Vk + 1(5; t+i + VkSik-i 



+ msik. 



where the function (pi^y) is introduced in Eq. (IB.lbj) . Notice that there is 



no 



second order contribution in the driving to the matrix element yll\t). This 
observation will be important later on. We now look at the matrix elements 
of y on Floquet states of the driven nonlinear oscillator from App I. We define 
for the following: 



yikit) := (0Ki^)|y|0fcW) = I]exp( 



in) 

Ik y 



(38) 
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where 



(n) 

Vik 



I^ex JO 



exp{+inujext){4>i{t)\y\(f)k{t)) = {{<pi,n\y\<Pk,G))- (39) 



We also define: 

yik,mn{t) ■■= {(pi,rn{t)\y\(t)k,ni^)) =exp{-iuextin-m))yikit). (40) 

We calculate the matrix elements within App I, which is exact in the driving 
amplitude. In particular from Eq. ([8]) we obtain: 

exp(—2iLUext)icQo + Cii) + exp(2zci;ea;^)(coo^ + c[i \ 
exp{-iujext), 

e(t) + ^ [2cOs{Uext){c^i^ + cSo'^) +2cOs{3Uext){c^ill + 4'^) 
App I 



2/oToi (t) = ^ 



ytr'nit) 



yZlit) 

e(t) + ^ [2cos(a;ext)(4l^ + 4'^) + 2 cos(3a;ext)(4? + 4'^) 
V2 L 



+2^2 cos(a;ex t)(4i^ + 4i^^) + 2y2cos(3a;ext)(4? + 



(3) , ^-3)^ 
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In the last derivations we used the coefficients c^^ introduced in Eq. ([13 
and the symmetry relations: 



In) 
^jk 



^kk ' 

_I-n) 
Hj ■ 



3^k 



(41) 



Inserting the actual form of the coefficients cj^^ we obtain: 



App I 

l/oo,oo 



App I 

?/io,io 



App I 

yii,ii 



COs{uJext) 



1 + 



3Fayl 



2M{u;l-n^), 



V2 



8hn 

COs{Uext) 



1 + 



9Fayl 



2M{ujl-n^)^ 



(42) 

(43) 
(44) 
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We observe that, as it is well known from the driven linear oscillator, terms 
of order or higher are absent at zero nonlinearity. The nonlinearity in- 
troduces a correction of which we neglect in the following. This is 
the reason why we neglected in Eq. f p2|) the second order contribution in 
the back transformation, as we can not rely on the F"^ contributions of the 
eigenstates originating from App II (see discussion in section H] concerning 
the limitations of App II and Appendix E), as long as these do not coincide 
with App I. 



5.2.2. yik{t) in the perturbative approach App II 

We calculate yik(t) for the one-photon resonance for the case j = 0. In 
the following we define |0o,o(^)) = (^|— 0,0)) and |0i,i(t)) = (t|+o,i))- Then 



00 



^exp(-inuJext){{-o,n\y\-o,o))- 



(45) 



To proceed we use Eqs. f p3|) and fl34j) for j = 0. Using the relation y = 
-^(a + a^), where a and a"'' are the annihilation and creation operators of the 
linear oscillator, we calculate the matrix elements: 



((0, n|(l + i?V»t(a + at)(l + RVf)\0, 0)) 



ypF nf(0) 

272^0,0-^1,-1 

A__(F)(5„„i + 5„,_i), 



((l,n + 1|(1 + RVpYia + at)(i + RVf)\\1, 1) 
n?(l) , nj{l) 



+ — + 

-E'1,1 ~" -E'2,2 -£^1,1 — -^2,0 



ypF 
2V2 



n?(0) 



(5n,l + 5n,-l)) = A+4F){6n,l + 6n,-l), 



{{l,n + 1|(1 + RVF)\a + at)(l + i?V»|0, 0)) 
((0, n\{l + RVF)\a + a)){l + RVf)\1, 1)) 



ni(0)5„_i = v4+_5„_i, 
^i(0)5„,i = v4_+5„,i. 

(46) 
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Note that = = 7^1(0) is independent of the driving. Consequently 
we find the result: 

l/oo!oo"(^) = Vot^^ exp{-iu}ext) +y!m^^ exp{iu}ext), (47) 

with 

»Sr' = !/So" (48) 

The other matrix elements are obtained in the same way. We give only the 
results: 

ytiji'\t) = = ^exp(-ma;ea:t)((+o,n+i|j/|+o,i))(49) 

n 

with 

sit" = y'n' (50) 
^ ^ (=in^ (F) + sin I COS + cos^ (F)) , 

Z/fo!ro^^W = (01,lWI^I<^O,oW) = X]^^P(~^^^ex^)((+0,n+l|y|-0,o))(51) 

= y[o^^ exp{-iujext) + 1/io exp(i(x;ea;t) 



with 

yir* ^ ^(-sin^|A_++sm|cos|[^__(F)-^++(F)l),(52) 
yfe" ^ -|(+si„|cos|!l^„OT->l.+(f)l + oos^|A,„), 

and 

n 

= yjt^^ exp{-iu}ext) + yoi^^ exp{iu;ext) 
17 



with 



y[o'\ 
-yit''- 



(54) 



Hence, yto,io^\t) = (yoi,oi^\t)^ • 

To compare we now expand the matrix elements in the driving strength and 
in the nonhnearity up to first order, using: 



sm 



cos 



Hi 

2 
Vj 



yoFni{j) 



1 yiF'nlU) 



1 



2 ■ 16 (^,+i,+i - Ej^oy 
Consequently Eqs. (1471) - ( l54l) yield: 

F 



+ 0{F\a'), 
+ 0{F\a') 



(55) 
(56) 



yZo\t) 
ytrii'\t) 



+0{a\F'), 

yo A 3ayo\ 



1 + 



3Fay^ 



tUeJ) + 0{c?, F^ 



+0{a\F^), 



1 + 



9Fayl 



2M{uj\ 



2 - D? 

ex 



(57) 

(58) 
(59) 



Hence App I reproduces the expressions for the matrix elements yik{t) ob- 
tained in App II up to first order in the driving F near the one-photon reso- 
nance. Note that, as in App I, the difference between yQQQQ^{t) and yiiii^{t) 
in the nonlinear contribution arises due to the contribution of states neigh- 
bouring the (quasi)-degenerate states. App I is not valid at resonance, as 
the generated degeneracy is not included in the perturbative treatment. For 
App II the resonance condition is essential for generating a doublet of degen- 
erate levels, requiring a finite nonlinearity. However the resulting energies 
and matrix elements can be expanded in the off-resonant regime and coincide 
with the results obtained within App I. Combining both approaches we can 
describe the whole range of possible driving frequencies, using Eqs. fH7|) -f l5^ . 
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6. Dissipative dynamics 



To include dissipative effects we use the system-bath approach introduced 
by Caldeira and Leggett |37|, where the bath is composed by harmonic 
oscillators with coordinate Xj momentum pj and frequency Uj and is cou- 
pled bilinearly to the system degrees of freedom. The interaction between 
system and bath is encapsulated in the coupling constants Cj. We have also 
introduced a counter term, proportional to y"^, which accounts for the renor- 
malization of the potential of the system due to the coupling. 
The total Hamiltonian of the system plus bath is given by: 



(60) 



where 



DO+B 



1 



— - — — nijUJ.x- 
2m,- 2 ^ ^ ^ 



1 2mja;| 



In order to specifiy the effect of the bath it is convenient to introduce the 
spectral density of the bath: 



J{uj) 



TC 



c2 



2 — ' rrijUjj 

7 = 1 J J 



(61) 



For a large number of bath oscillators the spectral density will become a 
smooth function and in the following calculations we use the Ohmic case: 
J{u}) = M'ju, leading to memoryless friction and white noise in the classical 
limit. We will assume for further calculation that the Duffing oscillator (DO) 
and bath are uncorrelated at time t = 0: 



Pdo+b(0) = Pdo(O) ® 



exp{-HB/kBT) 
tTB exp{-HB/kBT) 



(62) 



where Pdo(O) is the density operator of the Duffing oscillator at time t = 0. 
Because the bath consists of infinite degrees of freedom we assume the effects 
of the interaction with the DO system on the bath to dissipate away quickly, 
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such that the bath remains in thermal equihbrium for all times t. 
We wish to obtain an equation of motion for the reduced density operator 
pj3o(t) = trBPDO+B(t). Following EH, H Hsl a Born-Markov ap- 

proximation is applied and a Floquet-Markov master equation for the reduced 
density operator expressed in the Floquet basis of the DO is derived: 



(63) 



+ 



^ ^ exp[-i(n + n)Uext] [{N^a'-n + Npii',n')yal'y^p^pPa'P' 
«'/3' nn' 

T\T ('"■') ('"■) T\T ("■') (") 1 

where pap(t) = {(f)a(t)\pBo(t)\(pi3(t)) ■ A Lamb-shift contribution was disre- 
garded. The other quantities entering Eq. (163!) are 

Na(3,n = N{ea-e(s + hnuex), (64) 

^(lel), 



Nie) 



h 



-K(|e|) + 0(-e)], 



where 6{x) is the Heaviside function and ritt is the bosonic thermal occupa- 
tion number rithU) = 



coth 



2kBT 



. Additionally, y^^ are the Fourier 



coefficients defined of the matrix elements calculated in Sec. 15.2.11 , see Eq. 
(!39|) . For additional simplification we perform a moderate rotating- wave ap- 
proximation consisting in averaging the time-dependent terms in the bath 
part over the driving period T^^^ = 2tx jiOex [i^, lisf : 



Pafl(^) 



aP^a'P'Pa'p'-: 



(65) 



where p indicates the time average and the dissipative transition rates are: 

^a/3,a'/3' = ^(iVaa' -n + A''/?/?' -n)l/Sl/^7^"^ (66) 



The reason for its application is that it yields a time-independent stationary 
solution 



limf_s.oo Pa/3(t), which can be obtained by solving the linear 



system of equations: 

= --(Cfv - (^p)p% + ^ £a/3,a'/3'P^*/3'- 



(67) 
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7. Observable for the nonlinear response 

Using Eq. (l67j) corresponds to describe the long time limit, where the 
system has already reached the steady state. We calculate the mean value 
of the position operator in the stationary state, the nonlinear response: 

(^),t = limtr{^DoW^} = J]Pa>/3a(t), (68) 

a/3 

where tr{} denotes the trace over the oscillator degrees of freedom. Upon 
focussing on driving frequencies near and at the first resonance, the solution 
in the long time limit for a driven system is given by: 

{y)st ^ A cos{uext + 0), (69) 
where higher harmonics have been neglected. We introduced the amplitude: 

A = 2|5^p^,y(:^)|, (70) 

a/3 

and phase shift 

ImEa/3Pa>£'^ 



+ne I -ReJ^Pa^y^^a I - arctan 

al3 



Z^a/3 PapV^Sa ^ 



(71) 



7.1. One-photon resonance using PSA 

We solve the master equation close to resonance in the low temperature 
regime, imposing a partial secular approximation (PSA) which amounts to 
take only coherences of the involved resonant levels into account. For an in- 
termediate damping strength, i.e. the damping is of the order of the splitting 
of the one-photon resonance or larger, the system of equations is: 

= £oo,ooPoo(i) + 'Coo,iiPii(t) + 2£oo,oiRe [poi(i)] , (72) 

,llPll + -^01,01^01 + -^Ol.loPoi' 

Pii = 1 - Poo, 
where for simplicity we omit the superscript 'st'. 

If the damping 7 is smaller than the smallest splitting we can even neglect 
the coherences completely: 

= y~^>^aa,/3/3P/3/3- (73) 
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Note that we use as shorthand abbreviations: eo := eg with eigenstate |0) = 
|0o,o(^)) and ei := Cq with |1) = |0i respectively. Solving the system of 
equation for the stationary solution, we obtain: 



Poo — {~'^oo,ii['^oi,oi ~ -^01, 10 + (^0 — ^lY] + 2£oo,oi'^^oi,ii('^^ 

~'Coi,lo)} / {('^00,00 ~ '^00,ll)['Coi,01 ~ -^oi.io + (^0 ~ Cl) 

— 2£oo,oi('^oi,oo — 'Coi,ii)('^^oi,oi — -^oi.io)} , 

(-^oi.oi — '^^oi,io)['^^oi,ii + (-^oi.oo ~ -^01,11)^00] 
2 

01,01 



01,01 
21 



f^epoi - 1^ —r 

-^01,01 ~ ''-OKio + 1^0 - eij 



Impoi = jir-^ 77 rRepoi- (74) 

l''-01,01 ~ ''-oi.ioj 

To simplify the rates and obtain analytic results we restrict to low tempera- 
tures, ksT <^ huex- Moreover, we consider the vicinity of a resonance, such 
that \ea — C/^I is of the order of the minimal splitting, proportional to F, 
and within the limit of validity of App II, we obtain \ea — c/sl ^ ^ex- 
This corresponds to consider only spontaneous emission (see Eq. (IM|) ): 
Nai3,^i - J{\uJex\) and N^fi^i ~ 0. 
So we obtain: 



with J(\bJex\) = M'y\Uex\- 

7.2. One-photon resonance versus antiresonance 

We use these approximate low temperature rates and solve the master 
equation near the one-photon resonance to calculate the amplitude and phase 
of the steady state. We calculate the amplitude A in Eq. (170]) for the one- 
photon resonance using the energies Eq. (1281) up to second order in the 
driving and first in the nonlinearity and the expectation values Eqs. f H7|) - 
flM|) up to both first order in the driving and in the nonlinearity. 
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1.2.1. The amplitude in lowest order of the damping 

First we start with the analytical result for very low damping, where we 
only have to solve Eq. fl73|) resulting in: 



"'^00,11 2/01,1 , N 

Poo = 77 ^ r = ( 2 ^ 2 ^ + ^(7) (76) 

I't-ocoo — ''-oo,iiJ 12/10,1 + 2/01, iJ 

Pii = l-poo + C'(7)- 

Therefore the amplitude Eq. (!70!) in lowest order, i.e. zeroth order, in the 
damping is determined to be: 



A 



2 I [yo (2 sin r]oAA^^ cos r^o^^ _ (F) (77) 
+ ((5cos(277o) + 3)Al_ - AA+4F) sm{2r]o)A+^) A__(F) 



+6Al^A++{F) sin^r^o -2A 



sm 



I 2V2 ((cos(2r/o) + ^)Al_ + 2(A__(F) - A++(F)) sin{2r]Q)A^ 
+2(A__(F)-A++(F))2sin2(r/o))]|. 

The actual form for the amplitude is shown in the Figure [2l Interestingly, 
an antiresonance occurs, as already predicted in 22 



For comparison we show the result obtained by Peano et al. 221]. We 



see that both curves show an exact antiresonance, but the dip position is 
slightly different and our result shows an asymmetric line shape. The reason 
for the differences are explained in the following. The condition for an exact 
antiresonance, A 



n • (1) 

0, IS Poo 2/00 



(1) 

-P112/11 • 



As the result of 22 



IS given 



in lowest order without back transformation, the dip position occurs at the 
one-photon resonance, when the driving is such that sin ^ = cos ^ = 1/ \/2. 



This corresponds to resonance in 22| due to poo = Pu = \ and with 2/00 
—y^i = — sin Y cos Y"^' yielding a symmetric shape of the amplitude, seen 
in Figure [2l When all contributions linear in a are included, we find at 
resonance: 

(A+_+A__(F)-A++(F))2 , 1 



Poo 



2[Al_ + {A__{F)-A^^{F)Y]^ 2' 



(78) 



and 2/00^^ 7^ ~2/ii^- Moreover 74^y=^/4 7^ 0. The reason for not obtaining 
an equal population of the involved levels at resonance is due to the back 
transformation leading to the dressing of the eigenstates by states outside 
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Figure 2: Amplitude A at low damping (7 ^ yoFni{0) / \/2) versus the driving frequency 
Uex- Chosen parameters are: fl = 1, y^a/{hn) = 4/3 • 10"^, yoF/{hn) = 2^/2 . IQ-^, 
ki,T/{hn) — 10~^ and 7/r2 = 10~^. For comparison we give also the result obtained by 
Peano et al. in Eq. (43) in p^ . 



the quasi-degenerate doublet. In general, the dip position is determined by 
minimizing the amplitude A, Eq. f l77p . with respect to the external driving 

frequency. As the driving enters in both A and in t]q, the amplitude 

acquires a nontrivial cjex-dependence, such that the minimization can only 
be done numerically. The antiresonance does not occur exactly at rjQ = n /4, 
but it is shifted to little larger values (deviation oc 10*^^''-') with respect to 
the resonant case. This deviation arises due to the interplay of terms involv- 
ing trigonometric functions of rjo and explicitly driving-dependent prefactors, 
y4++(F) and A {F), resulting from the back transformation. 

7.2.2. Solution for the amplitude including higher orders in the damping 

We compute the amplitude in the low temperature regime, for fixed driv- 
ing amplitude and varying damping strengths by solving Eq. ( 17^ for the 
one-photon resonance Eq^ = Ei^i . Depending on the actual value of the 
damping, introduced by the bath, either an antiresonant behaviour, at small 
damping, or a resonant one, at large damping, is observed. The resonant/ an- 
tiresonant lineshape depends on the ratio of damping and minimal splitting: 
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%\/2'^j [l — g^a] -Fyo- In case of high damping, the amphtude is very small 
and broad showing resonant behaviour, where the asymmetry is smeared out 
completly. If we decrease the damping, the amplitude increases and the width 
shrinks until we reach a critical value for the damping. This critical value 
occurs when damping and minimal splitting are almost equal. Lowering the 
damping below the critical value, a cusp-like profile arises: the antiresonance. 
The treatment of higher orders in the driving and the interplay of driving 
and nonlinear ity introduce additionally an asymmetry in the response with 
respect to uj^x = + J^«l/o- ^^^y large damping we observe the cor- 
responding linear response (LR) of a linear oscillator with eigenfrequency 
^ + J^ct|/o) indicated by the black dashed line in Figure [21 So we can make 
the connection to the linear response of a driven damped harmonic oscillator, 
which is resonant at the frequency uj^x = + J^^l/o' if considering first order 
perturbation theory in the nonlinearity. 
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Figure 3: Amplitude A versus driving frequency at low temperature ksT/ (hil) = 10 ^ 
for different damping strengths 7. For the rest of the parameters we take yQa/{hil) = 
4/3 • 10-3, yoF/{hn) = 2^/2 . 10-5, n = l and varying damping: -f/n = 10"^ 5 • 10"^ 1 • 

lo-^ 3 • lo-^ 5 • lo-^ i • lo--^, 2 • 10-*. 



7.3. Phase 

The phase for the one photon resonance is given in Figure HI It jumps 
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1.00090 1.00095 1.00100 1.00105 1.00110 

a>ex (units of fl) 

Figure 4: Phase (j) versus driving fequency at low temperature kBT/{hVl) = 10^^ for 
different damping strengths 7. For the rest of the parameters we take a = 4/3 • 10"'^, 
F = 2^/2 . 10"^ = 1 and varying damping: 7 lO""^, 5 • 10^^ 1 • 10"^3 • 10"^ 5 • 
10-5,1 • 10-^,2-10-'^. 

by vr at the one-photon resonance and the step becomes smoother the larger 
the damping is. 

8. Conclusions 

In this work we have discussed two perturbative approaches, App 1 and 
App 11, to calculate the Floquet quasienergies and states of the quantum 
Duffing oscillator beyond a RWA approach. Additionally, the stationary 
dynamics of the expectation value of the oscillator position was obtained. 
Specifically, analytical results where derived off and near resonance with 
App I and at the one-photon resonance within App II. For App 11, based 
on Van Vleck perturbation theory, we also assumed that the driving is much 
weaker than the nonlinearity. We showed that in the parameter regime near 
resonance the Van Vleck approach recovers the results of App 1 based on the 
exact Floquet states of the driven linear oscillator, for both the eigenfrequen- 
cies in second order in the driving strength F and the matrix elements of the 
position operator in the Floquet basis to ffist order in F. The comparison 
allows to treat the quantum Duffing oscillator for arbitrary frequencies, as 
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well as an estimation of how good the Van Vleck approach is as App I is 
exact to all orders in the driving. As an application of our formalism we con- 
sidered the dynamics of a quantum Duffing oscillator coupled to an Ohmic 
bath and calculated its response near an one-photon resonance. Dissipative 
effects strongly affect the behaviour of the Duffing oscillator in the resonance 
region. We observe upon variation of the damping strength a transition from 
antiresonant to resonant behaviour and that the shape of the antiresonance 
with respect to the one-photon resonance condition is asymmetric. 
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Appendix A. Floquet theory 



The Floquet theorem states that a Schrodinger equation involving a time- 
periodic Hamiltonian {H{t) = H{t + T^^^) with period T^^^ = 27i/uji,x) (ssl . 



26|, l27|, |28|, |43 



zhdt\m) = H{t)m)) 

has a complete set of solutions: 

|^,(t))=exp(-ze,t/;i)|0,(t)), 



(A.l) 



(A.2) 



where \4>j(t)) 



,(t + T^^^)). The quasienergies and Floquet states \4>j{t)) 



solve the eigenvalue equation of the Floquet Hamiltonian: 

d' 



mm)) 



Hit) - ih 



dt 



(A.3) 



Defining \(pj,n{t)) = exp{—inUext)\(pjit)) and inserting into Eq. flA.3l) we find 
that \(j)j^n{t)) is also an eigenstate of the Floquet Hamiltonian, but with the 
eigenvalue e 



Sj — nhuex differing by a multiple integer of hwex- These 
Floquet states are therefore physically identical. In other words the spectrum 
of the Floquet Hamiltonian has a Brillouin zone structure, each Brillouin zone 



27 



being of size ujex- To find a complete set of solutions {\il^j{t))}, it is sufficient 
to consider only those Floquet states and quasienergics which lie within a 
single Brillouin zone, i.e. —fvjJex/'^ < < Moreover ej = ej^ and 

|0j(t)) := \(})jfi{t)). The eigenstates lipjit)) as well as the Floquet states 
are elements of the Hilbert space TZ. For convenience we also define 
T as the space of T^^^^ -periodic functions with the inner product: 

{f:9) = TfT- / dtnt)g{t). (A.4) 

-'■LOex Jo 

An ort honor malized basis of T is given by the functions: ipn{t) = exp{—inoJext), 
n integer. The basis set {(fin} is orthonormalized and complete: 

i^n,Vn) = Kw, (A.5) 

-E^nW^n(t') = 6T^Jt-t% (A.6) 



T 



where 5t^^^ is the T^^^^-periodic delta function. The scalar product in the 
composite Hilbert space TZ^T is then given by: 

((0,10,)) := -- / dt{<l>j{t)\Mt))- (A.7) 

The decomposition of \(pj(t)) into basis functions (fin{t) is equivalent to an 
expansion in Fourier series: 

Ihnit)) = J]exp(-i/a;e.i)|(/.5-^-"^), (A.8) 

Wn = 7fr- / dteM^nu,xt)\'PM (A.9) 

-l-Utex Jo 



2 fTujex 

—- / dt\4>j,-n{t)). 

-'■Wex Jo 



For a basis independent notation we introduce the state vectors \n) with 
(t\n) (Pnit)- Then, in the composite Hilbert space TZ <^ T we define the 

l</'.,n)) = El'^r"^)®IO' (A-10) 

I 
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and \(pj,n(t)) = (t\(j)j^n))- In particular, due to the orthogonality of the Floquet 
states \(f)j^rn{t)), it also holds the orthonormality relation: 

{{4>i,n\4>j,m)) = Tf. / dt{(f)i^n{t)\(l)j,m{t)) = Sij6nm- (A.H) 

-l-UJex Jo 

Additionally we define: 



^ {{<P^An<pJ,m)) = Tf^ dt(0,„(t)|H(t)|0,,„(t))-(A.12) 

Hence through the expansion of the Hilbert space it is possible to treat the 
time-dependent problem (lA.ip as a time-independent one by expressing the 
Floquet Hamiltonian in a basis of 7^ ® T. Which basis of 7?. T is the most 
convenient to express "H depends on the specific problem. For example in 
the basis {|0j>))} of the Floquet states the eigenvalue equation reads: 

ii\(j)j,m)) = (^j,m\4>j,m))- (A. 13) 



Equivalently, using the expansion f lA.lOp it also follows from f lA.12p and 

(EUD 

J2 {h^'-'^ - I'nuJeJu) |0f -"^^) = e.-.n.l^r'"^), (A.14) 
v 

where H'^^^ are the Fourier components in the Fourier expansion of H{t). 
The eigenenergies and eigenstates of l-L are known only in very few cases, 



among which the case of the driven linear oscillator, see [Appendix B[ For 
a generic time-periodic Hamiltonian only approximated solutions of Eqs. 
flA.12p or (IA.13P can be found. Two complementary approximation schemes 



are discussed for the Duffing oscillator case in section [2J 

Appendix B. Floquet states of the driven harmonic oscillator 



Following Husimi, Perelomov and Breuer [30|, |31|, |4J] the quasienergy 
spectrum can be determined exactly for a periodically driven harmonic os- 
cillator. The solution {y\ipj{t)) = ipj{y,t) of (]A.3|) for the time- dependent 
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Hamiltonian Hi^o{t), see Eq. (12a|) . is: 



i^jiy^ t) = <Pj{y - ^{t)) exp <^ -i 



1 1 



My) = (y\:j)o 



ex -— ] H- ( — 



(B.la) 
(B.lb) 



with |j)o the eigenstates of the hnear, undriven oscillator and Hj the 

Hermite polynomials. Finally S,{t) is the steady state solution of the corre- 
sponding classical equation, 



um+un^m = fit), 

which is for a driving of the form f{t) = —Fcos{uJext)'- 



(B.2) 



(B.3) 



The quasienergy spectrum for the harmonic oscillator is for a cosine-like 
driving term: 



,(0) _ M 



nnij 



+ 



2 J 4M{ul-n 



2\ ■ 



Finally, 



<pf{y,t) = 0^.(1/ -e(t))exp(-^^(F^^)) 



with 



1 

2h 



dt'fimt') 



x. 



i^ex J 



dt'f{t')m 
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Appendix C. Fourier components 

In this appendix we evaluate the matrix elements of the perturbation 
on the Floquet basis of the driven linear oscillator as well as the Fourier 
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coefficients 



a 



o{Mt)\Vam))o = J I dy'dy\{Mt)\y'){y'\f\yl{y''\<f>jmc-'^) 



a 



-oo 

CXD 



dyMy) [y + m]' Hv) 



oo 

oo 



poo 

7 / dyMy) [y' + ^my' + ^^'y 

J —oo 

+m'y+m'Wy) 

^ exp(-inujext) p(0fc(^)|K|0j(O)o"/ 



(C.2) 



(n) 



The Fourier coefficients are given, using the notation = -o."^) ' 



by: 



(0) 



4 



4, Q(2j + l)2/o'A^ + + 1) + \)yt + (C.3) 



+4.+4f va+i)a+2)(i+3)(i+4) 
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,,(±2) 
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+4,,+2^yo4V(j + l)(j + 2) 
+6kj-2^ylAyj{j -1)], 



(C.5) 
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Appendix D. Van Vleck perturbation theory 

In the foUowing we give a basic introduction of the Van Vleck perturba- 
tion theory [33|, |M, |35|, |36|]- It allows to calculate eigenenergies and eigen- 



functions of Hamiltonians H whose spectrum is splitted into well-defined 



manifolds (denoted by Greek indices) [34, |35|, l36|, |45|, |46|. Within Van Vleck 
perturbation theory an effective Hamiltonian H^s = exp{iS)H exp{—iS) is 
constructed whose spectrum is the same as that of the original Hamiltonian 

H = H^°^ + V, (D.l) 

but only connects almost degenerate levels within a given manifold. The 
eigenstates are then calculated from the eigenstates of the effective Hamilto- 
nian via a back transformation. For time-dependent Hamiltonians one intro- 
duces an effective Floquet Hamiltonian T^efr in the composite Hilbert space 
TZ®T. The effective quasienergies up to second order in the perturbation 
are determined by: 



((z, a|?iefr| j, o)) = Ei^a^ij + ((i, a\V\j, a)) 
+ \ 5^((*,«|V|fc,7))((A:,7im«)) ^ ^ 



(D.2) 



E., 



E, 



+ 



fc,7 



E. 



E. 



fc,7 



= E^^J,^ + {{i,a\V\3,a)) + {{i,a\{RV)^V\3,a)) 
+ {{i,a\V{RV)\3,a)), 
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where we introduced the reduced resolvent: 



R = Yl '\k,j)){{k,^\/{E-E,,, 

fc,7 



(D.3) 



The prime over the sum denotes that all the states belonging to the manifold 
a under consideration are excluded from the sum. The states of the original 
Floquet Hamiltonian are given by applying a back transformation connect- 
ing different manifolds: exp{-iS)\j, a)) = (1 - tS^^^ - iS^^^ + iS^^HS^^^ /2 + 
. . where: 



{{t,a\V\j,(3)) 



Ein — E 



(D.4) 



((z,a|\/|fc,7))((A;,7l^|j,/3)) 



E, 



E, 



A;, 7 



E, 



E, 



A;, 7 



E. 



k 



{{,.a\V\k,l3)}{{k,l3\V\j,m 



Ei R — Ei , 
1 



Ek,l3 — Ei^a 



{{z,a\V\k,a)){{k,a\V\j,f3)) 



The construction using the reduced resolvent is more easily comparable to 
conventional degenerate perturbation theory 47], |48||. After identification 
of the degenerate levels the modifications to the eigenvectors are given by 
calculating all possible matrix elements from the degenerate levels out of the 
manifold. 



Appendix E. Comparison for the states 

As seen in section 15.11 App II yields away from resonance the same 
quasienergies as App I expanded up to second order in the driving amplitude 
F. In the following we determine whether the states behave in the same 
way. For simplicity we compare the two approaches for the case of the driven 
linear oscillator. This corresponds in the perturbative approach, App II, to 
expand to zeroth order in the nonlinearity and first order in the driving. 

Appendix E.l. Floquet states in App I 

According to equation fIB.lap the Floquet states of the driven linear oscil- 
lator {|0j(t))o} can be obtained from those of the undriven linear oscillator 
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{|j)o} by applying a time-dependent translation: 



Um) = exp(^-^(a-at)J, (E.l) 

and accounting for a driving dependent phase exp{—i9{F'^ , t)), see Eq. (IB.6I) . 
Here d^ and a are the usual creation/annihilation operator related to the 
linear oscillator. To compare to the results of App II we give the Floquet 
states from App I up to first order in the driving amplitude and zeroth in 
the nonlinearity 

^7(e(^))expHe(F^t)) = l + ii^(dt-a) + 0(F2), (E.2) 

1/0 V 2 

\<PMo = Um)^M-^OiF',t))\j)o (E.3) 
= |j)o + 1^ + + l)o -Vj\J- 1)o) + 0{F'). 

In the composite Hilbert space 1Z®T the Floquet state corresponding to ej°Q 
is: 



+ VTTIb- + 1. +i))o - VIIj - 1, -i»o - VJb - 1, +i))o) ■ 

Appendix E.2. Floquet states in App II 

In App II the counterpart of Eq. (IE.4p is the state: 



f,o))o = |j,0))o + 



v(7TT|j + 1, +l))o VJTTIj + 1, -l))c 



v/j|j-l,-l))o , Vj|j-1,+1))o~ 



Note that by reducing to the linear case not only the energies Ej^ and 

but also Ej_i _i are equal if ^ Uex- These are the dominant contributions, 

as — OJex) » 1 if ~ UJex- 



,-,o))o=|j,0))o + 



v^7TT|j + l,+l))o v^|j-l,-l))c 



2V2MyQn \ UJex - ^ex-^ 

(E.5) 
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Appendix E.2.1. Comparison 

By comparing the states we see that the exact hnear oscillator Floquet 
states from App I are proportional to ^ (t) and hence are obtained as a linear 
combination of rotating and anti-rotating contributions. This accounts for 
the fact that the spectrum of the underlying oscillator is equidistant. Con- 
sequently the driving can excite both upwards and downwards transitions 
|j)o |j + l)o) |j)o |j — 1)0 with the same weight. In contrast, App 
II in a vicinity of a A^-photon resonance includes rotating and antirotating 
terms but with different weights. In fact, by choosing a resonance condition 
certain transitions are preferred, namely those within the resonant levels and 
those lying closest to these, giving the dominant contributions (see Eq. lE.Sp . 
The actual form for the states in App II is determined by the structure of 
the quasienergy-spectrum, i.e., by the resonance condition determining the 
structure of the manifold. As the manifold structure of App II is crucially 
depending on a finite nonlinearity to obtain a doublet, a reduction to the 
linear case is not possible. It would destroy the ordering of the manifolds, as 
not only Ej^ and but also Ej^n,n, n > 1 and Ej_n~n, n > 1, become 

degenerate. The way of treating the Floquet Hamiltonian thus results in 
different symmetry properties and weights for the corresponding states. 
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